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Abstract A bright prominence associated with a coronal mass ejection (CME) 
was seen erupting from the Sun on 9 April 2008. This prominence was tracked by 
both the Solar Terrestrial Relations Observatory (STEREO) EUVI and CORl 
telescopes, and was seen to rotate about the line of sight as it erupted; therefore, 
the event has been nicknamed the "Cartwheel CME." The threads of the promi- 
nence in the core of the CME quite clearly indicate the structure of a weakly to 
moderately twisted flux rope throughout the field of view, up to heliocentric 
heights of 4 solar radii. Although the STEREO separation was 48°, it was 
possible to match some sharp features in the later part of the eruption as seen 
in the 304 A line in EUVI and in the Ha-sensitive bandpass of CORl by both 
STEREO Ahead and Behind. These features could then be traced out in three- 
dimensional space, and reprojected into a view in which the eruption is directed 
towards the observer. The reconstructed view shows that the alignment of the 
prominence to the vertical axis rotates as it rises up to a leading-edge height of 
~2.5 solar radii, and then remains approximately constant. The alignment at 2.5 
solar radii differs by about 115° from the original filament orientation inferred 
from Ha and EUV data, and the height profile of the rotation, obtained here 
for the first time, shows that two thirds of the total rotation is reached within 
~ 0.5 solar radii above the photosphere. These features are well reproduced by 
numerical simulations of an unstable moderately twisted fiux rope embedded in 
external fiux with a relatively strong shear field component. 
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1. Introduction 



The bandpass of the inner coronagraph (CORl) telescopes on the two Solar 
Terrestrial Relations Observatory (STEREO) spacecraft runs from 650-670 nm 
( [Howard et a/., 20081 ). This range was selected to include the hydrogen Ha fine 
at 656 nm, so as to be sensitive to erupting prominences associated with coronal 
mass ejections (CMEs). A particularly bright prominence eruption was observed 
by both CORl telescopes on 9 April 2008. Although the relative contributions 
of Ha emissions and Thomson-scattered light cannot be derived from the CORl 
observations alone, it is presumed that Ha is a significant contributor for the 
extremely bright parts of the prominence early in the eruption. This CME and 
prominence eruption was also observed by the STEREO Extreme Ultraviolet 
Imager (EUVI) telescopes ( Howard et a/., 2008 ) before entering the CORl field 
of view, as well as by TRACE at 1 71 A ( [Handy et al, 1999D, and by the Hinod e 
XRT telescope ( [Golub et al, 2007[ ) and EIS spectrometer (IC ulhane et al, 2007 ). 
ILandi et al.\ (2010) analyze this event, combining Hinode/XKT and EIS data, 
together with Solar and Heliospheric Observatory (SOHO)/UVCS/EIT/LASCO 
and STERE0/EUVI/C0R1/C0R2 images, to characterize the thermal proper- 
ties of the ejected plasma, and constrain the heating rate. The post-CME current 
sheet following this event is analyzed by [Savage et aL\ (|2010p and IKo et a/.1 
(|2010p . IPatsourakos and Vourlidasl (|201ip determine the geometrical parameters 
of the CME and post-CME current sheet in the field of view of the outer coron- 
agraph (C0R2) on STEREO. This event has become known as the "Cartwheel 
CME" because of the highly visible rotation around the line of sight seen during 
the initial stages of the eruption. Rotation about the vertical direction can also 
be observed through triangulation as the prominence material rises through the 
EUVI and CORl height ranges. The plane of sky motions which give this event 
its "Cartwheel" name is actually a combination of the rotation about the vertical 
axis together with a defiection of the prominence in latitude and longitude 
during the initial rise phase. A rotation around the axis of the prominence, 
known as "roll effect" ( [Martin, 2003[ ), may have contributed as well. Further 
rotation about the vertical direction is indicated by the CME orientation in the 
C0R2 data. This paper presents an exploration of the rotation about the vertical 
direction in the EUVI and CORl height ranges, using triangulation from the two 
STEREO viewpoints. The height profile of the rotation angle for the erupting 
prominence is obtained, which could not be done prior to the STEREO mission. 

The separation of the two STEREO spacecraft at this time was 48°. Co- 
identification of features in both views is quite difficult with such a large sep- 
aration. However, it was possible to match some sharp features in the later 
part of the eruption as seen by EUVI in the He ll line at 304 A. Features were 
also tracked in white light (with a possible Ha component) as seen by CORl, 
although not necessarily the same features seen at 304 A. 
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Figure 1. Ho; image of the filament as seen on the disk on 5 April 2008, 01:51 UT, before the 
first eruption on that day, and on 6 April, 02:44 UT (circled), reformed after the two eruptions 
on 5 April (from the Yunnan Astronomical Observatory in Kunming, China). 



The derived height profile of prominence rotation is compared with the ori- 
entation of the polarity inversion line of the radial field component (PIL) at 
various heights, to check whether the rotation is simply the result of alignment 
with the PIL. The radial field component is obtained from a potential- field 
source-surface (PFSS) extrapolation. Additionally, the rotation and rise profiles 
of the prominence, as well as the STEREO images, are compared with the 
corresponding data obtained in a series of CME simulations in a companion 



paper (Kliem, Torok, and Thompson, 1201 1( in the following Paper II), in order 
to explore a broader range of possible origins of the large total rotation ob- 
served. The comparisons suggest that the major part of the rotation was due to 
the presence of a shear field component of considerable strength in the source 
region ( [Isenberg and Forbes, 2QQ7[) and that the weak helical kin k instability of 
a moderately twisted fiux rope ( Torok, Kliem, and Titov, 2004 ) contributed as 
well. Field line plots of the best matching case are included here to demonstrate 
the correspondence with the observed overall shape and helical threads of the 
rising prominence. 



2. Data Analysis 
2.1. Source Region 

The eruption on 9 April 2008 occurred in the remnants of NOAA active region 
(AR) 10989, located about 23° behind the west limb as seen from Earth, and 
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Figure 2. STEREO Ahead EUVI 195 A images of AR 10989 before and after the second 
filament eruption on 5 April 2008, displaying an area of 640^^ on a side. The image on the left 
at 19:45:30 UT shows the full extent of the filament on the souteastern end at the onset of 
the eruption (up to the middle of the circle), and the image on the right at 21:05:30 UT (with 
the circle plotted at the same position) shows that a dimming has developed at this location. 
Therefore, the erupting fiux should be rooted in this area. Comparison with Figure ^ shows 
that only the upper branch of the filament is visible in Ha. Overplotted grid lines represent 
Carrington longitude and latitude. 




Figure 3. STEREO Ahead EUVI 171 A image of the reformed filament on 6 April 2008, 
15:01 UT. The area shown is 320^^ on a side. The circles mark the ends of the filament material 
visible at this wavelength. Their distance is 175 Mm. The main body of the filament as seen 
in this line is oriented at an inclination of ?=i26° to the east- west direction on the Sun. 




Figure 4. Magnetogram of AR 10989 on 31 March (left) and on 4 April 2008 (right) taken 
by SOHO/MDI. The overplotted bars indicate rough estimates of the distance between the 
center of each polarity (in a center-of-gravity sense), which is an important input parameter 
for the numerical modeling of the event in Paper II. 
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had an onset time near 09 UT. The filament had been visible nearly through- 
out the disk passage of the active region, which was in its decaying stage and 
spotless after 31 March. The filament erupted twice on 5 April 2008, reforming 
afterwards. These eruptions were visable on disk to EUVI on STEREO Ahead. 
Both eruptions appeared to propagate inclined from the radial direction toward 
higher southern latitudes, similar to the eruption on 9 April. The inclination 
made it difficult to discern rotational motions based on images from only a 
single viewpoint. However, a stereoscopic reconstruction of the second eruption 
revealed a considerable rotation of about 90° ( |Bi et a/., 2011D . Interestingly, the 
rotation was in the clockwise direction, opposite to that of the eruption on 
9 April. 

As with most filaments, it is difficult to discern the magnetic connections at 
its ends. However, both eruptions on 5 April produced a dimming and endpoint 
brightenings in the EUV dWang, Muglach, and Kliem, 2009D near the southeast- 
ern end of the filament, which were located in negative polarity, the polarity 
which dominated the southwestern side of the filament channel, so that the fila- 
ment was dextral according to the classification of lMartinl (1998). This is further 
supported by the right-bearing orientation of the filament barbs (see Figured]). 
Dextral filaments are embedded in field of left handed chirality, as indicated 
by the typical skew of the overlying arcade of coronal loops ( [Martin, 199 81). If 
the assumption of flux rope topology holds for the filament, then the filament 
itself is also threaded by left-handed field. Otherwise, it may have the opposite 
chirality (see IMartin and McAllisterl [T9971 [Ruzmaikin, Martin, and Hu| [20031 
[Muglach, Wang, and Kliem, 2009J. If a filament in dominantly left-handed he- 
lical field rotates upon eruption, the direction of rotation is generally found to 
be in the counter-clockwise direction, as expected from helicity conservation 
( [Rust and LaBonte, 2005} [Green et al, 2007[ ). 

Note that the dextral chirality of the filament is an exception to the usual 
hemispheric rule ( [Martin, Bilimoria, and Tracadas, 1994[ ). However, [Pevtsov, Balasubramaniam, and Rogers] 
(2003) show that 16-25% of filaments do not follow the hemispheric rule, and 
that the rule is weaker for filaments in active regions. 

On 5 April, the filament was oriented at an inclination of ^ 24° to the east- 
west direction on the Sun and had a length of ^ 135 Mm in the image shown 
in Figured] (top). However, the dimmings and endpoint brightenings formed in 
the eruptions later on the same day suggest that the true extent of the flux in 
the filament channel was larger by about one third, 175 Mm, see Figure [2] 

Due to foreshortening, the Ha images on the following day show the filament, 
its magnetic connections and barbs far less clearly, but it can be seen that the 
filament reformed in a similar location (Figured! bottom). The better perspec- 
tive of the Ahead images indicates that the main body of the reformed 
filament was nearly straight, oriented at a tilt angle of ~ 26° to the east- west 
direction, and extended across a similar length as on 5 April (Figure [3]). This is 
most clear if the 171 A and 304 A images are viewed in animated format. As 
far as the increasing foreshortening allows a judgment, the filament appeared to 
keep its orientation and shape through the subsequent days until a new eruption 
launched the Cartwheel CME. 

Since the numerical modeling in Paper II reveals a strong influence of the 
distance between the polarities in the source region on the rotation of the ejected 
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flux rope, we attempt to estimate this parameter. Figure |4] shows magnetograms 
of the region from the SOHO Michelson Doppler Imager (MDI) on 31 March and 
4 April 2008. The overplotted bars run through the middle of each polarity in a 
one-dimensional center-of-gravity sense, providing an indication of the distance 
between the main polarities. This value increases from ~ 40 Mm to a range of 
~ (40-75) Mm over the four-day time span. A range twice as wide may have been 
characteristic of the configuration after the similar time span to the eruption on 
9 April. It appears impossible to estimate a single relevant value. We will adopt a 
distance of 90 Mm as a base value in most simulations in Paper II, but also study 
the influence of its variation. This base value lies in the middle of the estimated 
range, so it is consistent with the fact that much of the rotation occurred in 
the course of the radial ascent above the south-east edge of the remnant active 
region, which is characterized by a considerable distance between the polarities 
(see Section [23]) . 

The magnetograms in Figure |4] show an overall displacement of the two polari- 
ties in the direction of the PIL. While the negative flux is more or less uniformly 
distributed along the fllament channel, the positive flux is more concentrated 
near the northern end. This results in a net shear fleld component pointing 
along the fllament channel in the southeastward direction. Coronal fleld lines 
that arch above the fllament thus possess a skew corresponding to the upper 
part of a left-handed helix, in agreement with the implication made above that 
the fllament was dextral. 

We have also estimated the height of the prominence from data taken by 
the EUV Imaging Telescope (FIT) ( ,Delaboudiniere et al.^ 1995J onboard SOHO 
on 8 April, when the prominence was close to the limb. The material seen at 
that time in the 304 A line extended up to a height of ~ 0.04 Rq above the 
photosphere, and the absorbing material seen at 171 and 195 A extended up 
to ^ 0.033 R(T). From EUVI 304 A data about an hour before the eruption, we 
estimate the height to lie in the range ^ 0.05-0.06 Rq^ so it is clear that the 
prominence experienced a slow rise before the eruption. 

Motions along the prominence could be seen in the 171 A channel of FUVI- 
Ahead on 9 April from about 08:20 UT onward. The flrst upward motions of 
prominence material along the CME path can be discerned by comparing the 
images at 08:48:30 and 08:51 UT. ILandi et aL \ (|2010p give a rather conservative 
estimate of the CMF start time, 09:10 UT, but the 171 A images, taken at 2.5 min 
cadence, show that most of the prominence body visible to EWl-Ahead was 
already moving by 08:53:30 UT, the start time quoted in [Savage et aL\ (|2010p . 
Since the event commenced in the absence of a strong perturbation (no signs 
of a signiflcant brightening, of a jet, or of perturbations resulting from nearby 
activity were seen), it must have developed from a small perturbation when the 
conflguration was near the boundary between stable and unstable states. In such 
cases the initial motion of the unstable flux is expected to behave exponentially, 
which requires some time to develop to a level that causes visible changes. Rapid 
changes in the vertical position were flrst seen between the EWl-Ahead images 
at 08:48:30 and 08:51 UT. Therefore, the actual start time should lie before 
08:51 UT. 

A helmet streamer, best seen in the C0K2-Ahead images, extended above 
AR 10989 in the plane-of-sky projection. In order to estimate the orientation of 
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Figure 5. Location of the heliospheric current sheet for Carrington Rotation 2068 from a 
potential-field source-surface extrapolation. The current sheet runs nearly east-west above the 
estimated CME direction of propagation as indicated by the vertical bar. The green and red 
regions represent open field regions of negative and positive polarity respectively. 

the heliospheric current sheet above the streamer, and whether it was magneti- 
cally connected with the PIL in the active region, we ran a PFSS extrapolation 
and a Wang-Sheeley-Arge model for Carrington Rotation 2068 at the Commu- 
nity Coordinated Modeling Center (CCMC), using a source surface location at 
2.5 solar radii for both models. Figure [5] shows the PIL and corresponding field 
lines originating at the 2.5 Rq source surface of the PFSS model, which yields 
a nearly east-west orientation of the heliospheric current sheet at the longitude 
of radial CME propagation estimated in Section 12.31 The Wang-Sheeley-Arge 
model gives the same result. It is clearly seen that the field at large heights 
above the active region is dominated by the large-scale field of the Sun, which 
points essentially opposite to the shear field component in the active region at 
small heights. 

A similar representation of the radial field component of the PFSS model at 
several representative heights is given in Figure [6l The PIL above AR 10989 
is part of a magnetic structure which is distinct from the PIL at the base of 
the heliospheric current sheet, and which is seen as a separate ring in the low 
corona {e.g. at 1.15 Rq). With increasing height, the ring-shaped PIL disappears 
(between 1.2 and 1.3 solar radii), and the PIL associated with the heliospheric 
current sheet gradually approaches the equator at the estimated longitude of the 
CME. 

2.2. 3D Reconstruction 

The SolarSoft routine scc_ineasure was used to measure the three-dimensional 
location of features in the STEREO EUVI and CORl images. The action of 
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Figure 6. Polarity inversion line of the radial field component in the PFSS model, (a) in the 
photosphere, (b) at the estimated initial height of the prominence, (c) at the largest height 
that shows an orientation of the PIL above AR 10989 relatively clearly, and (d) at 1.3 Rq 
where the PIL above AR 10989 is completely gone. Dots in panel (a) mark the end points 
of the prominence shown in Figure [S] and the vertical line in panel (d) marks the estimated 
longitude of the CME, as in Figure [S] 



this program is as follows: The user is presented with two side-by-side images, 
one from each of the two STEREO spacecraft. The images are selected so that 
they represent the same observation time, with a slight offset to account for 
the difference in light travel time from the Sun. The user can zoom in on the 
region of interest in the two images, and adjust the color table and data range 
to optimize the appearance of the feature being measured. A point is selected 
on one image with the cursor. The program calculates the three-dimensional 
line of sight represented by this point, and then overplots the projection of 
this line onto the image from the other satellite. This is known as an epipolar 
line. Since both the EUVI and CORl optics produce a gnomonic projection on 
the CCD detector, straight lines in space will always appear as straight Hues 
in the image. The feature selected by the user in the first image must appear 
along the epipolar line drawn by the program in the second image. The correct 
location along this line is selected by the user, which leads to another line of 
sight calculation which intersects the original line of sight. The intersection of 
these two lines determines the three-dimensional (3D) location of the feature. 
In this investigation, we use this technique to find the 3D locations of multiple 
points along prominence threads. 

The main difficulty in applying this technique to the prominence eruption is 
source confusion. On 9 April 2008, the two STEREO spacecraft were separated 
by 48° along the ecliptic plane. Thus, the appearance seen by STEREO Ahead 
was quite different from that seen by STEREO Behind. This made it very difficult 
to locate features which could be positively identified to be the same in both 
images. Another goal was to identify features which could be tracked through 
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several frames. (However, no attempt was made to match identifications between 
EUVI and CORl.) In spite of these difficulties, it was possible to identify several 
features which could be identified in both views, and which could be tracked 
through several frames. 



2.3. EUVI and COR Measurements 



Figures [7] and [8] show the tracked features as seen in the 304 A channel of EUVI- 
Ahead. The colored lines represent the filamentary features that were measured 
at each time step. Each line has a different color, and a dot at one end, to 
better show the relationship with the reprojected graphs along the right side. 
The different colors do not necessarily represent the same threads at each time 
step. The He ll 304 A channel was chosen as being most representative of the cool 
prominence material seen later in a combination of Ha and Thomson-scattered 
light by CORl. Although the eruption is clearly seen by Ahead at 9:06 UT, and 
by Behind at 9:36 UT, source confusion made it impossible to track features 
earlier than 10:06 UT. Initially, little sense could be made from an examination 
of these data viewed in 3D {e.g., through an anaglyph representation). However, 
because the measurements are made in 3D, they can be reprojected to another 
viewpoint. By selecting the proper viewpoint, it's possible to see aspects of the 
data that are not evident in the frame in which the data are taken. The plots 
to the right in Figures [71 and [8] show the same data as in the images, but from a 
viewpoint at Stonyhurst longtitude 98° west (relative to Earth), and 24° south. 
It is estimated that from this viewpoint, the CME is traveling straight toward 
the observer. 

From this perspective, a certain organization can be perceived. The primary 
alignment in the 10:06 UT data is vertical, i.e. aligned north-south, with the 
upper prominence threads lying almost directly over the lower threads in the 
frame of the figure. However, in the later images the primary alignment is 
distinctly sloped, with the upper threads lying toward the east, and the lower 
ones toward the west. This is particularly well seen in the image at 10:46 UT in 
Figure m where the data display an overall slope of —50°. As one goes through 
the time steps, a counter-clockwise rotation can be perceived as one steps from 
10:06 to 10:46 UT. 

The prominence structure must have rotated through ^65° from its original 
orientation to achieve the nearly north-south structure seen at 10:06 UT, followed 
by an additional ^ 50° to reach the maximum rotation at 10:46 UT. Thus, the 
full amount of rotation experienced by the prominence in the EUVI height range 
was on the order of 115°. 

The EUV observations of the early part of the eruption also show a shift of 
the entire prominence structure southward, from about -14° to -24° latitude. 
Compared with the earlier Ha observations (Figure [1]), the prominence must 
have also shifted ^ 15° eastward in the early phase of the eruption. These 
initial motions were followed by radial propagation at the new latitude and 
longitude, ^ 98W24S. This is consistent with the analysis of [Savage et a/.| (2010) 
and IPatsourakos and Vourlidas[ (|201ip . The latter authors fitted a croissant- 
shaped flux rope model to C0R2 images from both STEREO satellites and 
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Figure 7. Images and plots of the prominence eruption from 10:06 to 10:26 UT, as seen in 
the 304 A channel of the STEREO EUVI telescopes. Subsequent frames are shown in Figure[8l 
Along the left are shown the KUYl-Ahead images with the tracked features overplotted. The 
lower-right corner of each image is outside the EUVI field of view, causing the threads to 
appear truncated in height at 10:26 UT and subsequent times. The plots along the right show 
the same data reprojected to a viewpoint at Stonyhurst longitude 98° west (relative to Earth) 
and latitude 24° south. Axes are in units of solar radii. 



found that the radial propagation of the CME projects back to a position 11° ±5° 
behind the hmb, as seen from Earth, and 17° ±3° south. Savage et al. (|2QlQp 



quote a hehographic position of 113° for the associated active region at the 
time of the eruption, and describe the initial motion of the prominence as being 
toward the observer. The measurements shown in Figures [7] and [8] were derived 
from the radial phase of the eruption. 

The same analysis performed on EUVI was also performed on the CORl 
images. At the location of the eruption, the CO^l- Ahead occulter edge is at 
1.37 solar radii, and the CO^l-Behind edge is at 1.59 solar radii. The features 



SOLA: prominence_final.tex; 16 December 2011; 1:11; p. 10 




Figure 8. Continuation of Figure[7]for time steps 10:36 to 10:56 UT. 

seen by CORl are a continuation to higher radial distances of the truncated 
304 A features seen by EUVI. There is a smah region of overlap between EUVI 
and CORl on both spacecraft, and the observations agree within this region. 
The CORl results are shown in Figure [9l The prominence is seen to maintain 
the orientation seen in the final EUVI images. There is a slight indication of a 
shift back in the clockwise direction. However, the orientations for the final two 
time steps at 11:15 and 11:25 UT are based on only parts of the prominence 
legs, so the apparent reverse rotation may not truly represent the behavior of 
the loop as a whole. 

IPatsourakos and Vourlidasl (j2Qlip estimated the orientation of the CME flux 
rope from their croissant model fit at a leading-edge heliocentric distance of 
about 13 Rq. They found that the flux rope (croissant) axis was inclined to 
the east- west direction on the Sun by a tilt angle of —4° ± 7°. This indicates a 
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Figure 9. Images and plots of the prominence eruption from 10:45 to 11:25 UT, as seen 
in the CORl telescopes. Along the left are shown the COHl-Ahead images with the tracked 
features overplotted. The plots along the right show the same data reprojected to a viewpoint 
at Stonyhurst longitude 98° west (relative to Earth) and latitude 24° south. Axes are in units 
of solar radii. 
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further counter-clockwise rotation to a total value of 150° ± 7° from the original 
orientation estimated in Section 12.11 (The other possible interpretation of a 
clockwise rotation by ~ —145° following the counter-clockwise rotation in the 
EUVI height range appears far less likely.) It should be noted that the croissant 
model does not include any writhe of its axis, so that the uncertainties may be 
higher; however, this is not expected to change the result by a large amount. 
The rotation angle at large heights represents a close alignment of the CME 
flux rope with the heliospheric current sheet (Figure [5|), which corresponds to 
the suggestion in |Yurchyshyn| (|2008p and |Yurchyshyn, Abramenko, and Tripathi 
(120091) . 

To our knowledge, the rotation of the prominence by ~ 115° in the corona 
and by ^ 150° up to 13 Rq belongs to the largest ever inferred in this height 
range. [Yurchyshyn, Abramenko, and Tripathi| ()2009|) report one case of rotation 
by ~ 143° and six cases in the range 80°-100° out of a sample of 101 halo CMEs 
observed by SOHO/LASCO, i.e., up to a distance of 30 Rq. The largest rotation 
found at 1 AU is on the order of 160° ( |Dasso et a/., 2007| |Harra et al, 2007 ). 



It is also of interest to note that the CORl images resolve the basic structure 
of the dense prominence material in the core of the CME out to a heliocentric 
distance of four solar radii. The CME core has the structure of a weakly to mod- 
erately twisted flux rope: a single flux loop with threads that are systematically 
but only weakly to moderately inclined to the axis of the loop. These indications 
of twist prove to be the primary observational finding that sets a preference for 
a slightly kink-unstable flux rope above a kink-stable flux rope, which match 
the observed rise and rotation characteristics to a comparable degree in the 
numerical modeling of the event in Paper II. Only few observations have revealed 
a flux rope structure for an erupting prominence so clearly in this height range 
(see, e.g.. IPlunkett et al\ [2QQQ1) . The images also indicate a nearly self-similar 
evolution of the flux rope throughout the instrument's field of view of (1.5-4) Rq. 
The outer coronagraph C0R2 imaged the CME core as well but did not resolve 
the details of its structure. 



2.4. Rotation and Rise Profile 



Figure [TO] summarizes the rotation of the erupting prominence as a function of 
time from both EUVI and CORl. The values plotted were obtained by aver- 
aging the orientations of the individual threads whose 3D positions could be 
reconstructed. The values obtained at 10:16 and 10:26 UT are considered to be 
less reliable than the others in the figure, since at these times the threads which 
could be located in 3D do not appear to be organized in a well defined structure, 
as, for example, flux rope legs (see the right panels in Figure [7|). It was not 
possible to derive the 3D prominence orientation for times before 10:06 UT due 
to source confusion. However, the visible appearance of the prominence as seen 
from STEREO Ahead between 9-10 UT is consistent with considerable rotation 
throughout this time period, leading to the conclusion that most of the rotation 
occured during the initial phase of the eruption. 

We also used scc.measure to derive a time- height curve for the prominence 
eruption. This was done by matching points at the leading edge of the promi- 
nence in both images. Since these points in the two images are not guaranteed 
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Figure 10. Orientation of the prominence body as a function of time. The dashed hne rep- 
resents the original orientation as estimated from the FiUVl-Ahead 171 A image on 6 April 
in Figure [3] The points at 10:16 and 10:26 UT are considered to be less reliable than the 
others — see text. 

to represent exactly the same point in 3D space, the results can only be treated 
as an approximation. However, since the extent of the prominence material is 
relatively small, the approximation is quite valid. The resulting time-height curve 
for the top of the prominence is shown in Figure [Til The CORl values appear to 
be a bit higher than one would expect from an extrapolation of the EUVI data. 
That may represent a difference between the appearance of the prominence at 
304 A versus white light. In addition, the heights derived from EUVI data after 
9:56 UT may systematically fall short of the true heights. It is possible that the 
traceable threads did not extend up to the true top of the structure in these 
images, since an upward extension of the loop-shaped structure in the image at 
9:56 UT had faded by 10:06 UT. The estimated height at 10:26 UT is only a 
lower limit given by the edge of the EWl- Ahead field-of-view. 

A quadratic fit to the EUVI data gives an initial speed of 67 km s~^, accelerat- 
ing to 173 km s~^ as the prominence leaves the EUVI field-of-view. A quadratic 
fit to the heights measured by CORl gives a velocity starting at 206 km s~^, 
and accelerating to 379 km s~^ by the end of the CORl data. There's some 
evidence from the CORl quadratic fit that the acceleration strongly decreased 
when the leading edge of the prominence reached a heliocentric height of 3 solar 
radii. This is borne out by a linear fit to the CORl data points above 3 solar 
radii, which gives a slightly lower velocity of 327 ± 9.2 km s~^. The combined 
EUVI and CORl data in Figure [TTI give the visual impression that much of the 
acceleration actually occurred up to ^ 2.5 solar radii. The C0R2 data give a 
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Figure 11. Heliocentric heights of the prominence leading edge vs. time. The value for 
10:26 UT is shown as a lower limit, since the prominence is right at the edge of the EUVI- A/ieaci 
field-of-view at this time. 



speed close to 400 km s~^. All of these values refer to estimated positions of the 
apex point of the visible flux rope structure in the core of the CME. 

ILandi et al\ (2010) also derive a time-height profile for the leading edge of the 
prominence material, using a combination of STEREO Ahead and SOHO data. 
Their results are very close to those presented here, with the possible exception 
of the amount of acceleration. ILandi et a/. I quote an acceleration of 59.6 m s~^, 
while our results are more consistent with a smaller value of 37.9 ± 4.0 m s~^ 
over the time range covered by Figure [Til They also give the projected velocity 
of the CME leading edge in the plane of the sky of STEREO Ahead^ which 
approaches 700 km s~^, still rising, at a heliocentric height of 3 Rq. 

Combining the information in Figures [10] and [TTJ we obtain the rotation of 
the prominence as a function of its leading edge height, see Figure [T2j (As there 
are only three time steps in common between the EUVI data sets, we have 
interpolated the EUVI and CORl height data from Figure [11] to the times of 
the EUVI measurements in Figure [10] to better show the trend.) These data 
emphasize again that much of the rotation is acquired in the height range h < 
R(7)/2^ or ~2 initial footpoint separations, above the photosphere. 

2.5. Comparison With the Orientation of the PIL 

The comparison of the PFSS extrapolation with the orientation of the promi- 
nence at large distances from the Sun obtained in IPatsourakos and Vourlidas] 
(|201ip shows that the prominence has approached a close alignment with the 
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Figure 12. Prominence rotation vs. heliocentric height of its leading edge. The diamond 
symbols represent EUVI and CORl height data from Figure [TT] interpolated to the times of 
EUVI orientation measurements in Figure fTo] 



heliospheric current sheet (Section l2.3p . An obvious question is whether such 
ahgnment could also lead to the rotation in the coronal height range studied 
here. We will see in Paper II that the rotation of erupting flux ropes and the 
change of the PIL orientation with height always point in opposite directions 
if the eruption originates in a simple bipolar region. However, Figure [6] shows 
that the topology of the field above AR 10989 is more complex. In particular, 
the field component along the filament channel, which influences the rotation of 
erupting flux ( [Isenberg and Forbes, 2QQ7t , changes sign with increasing height. 
In Figure [6l it is clearly seen that the relevant section of the PIL keeps a nearly 
constant orientation up to about 1.15 Rq, above which it starts to shrink no- 
ticeably, disappearing completely by 1.3 Rq. Thus, any rotation in the height 
range up to about 1.15 Rq could not be caused by a changing PIL direction. In 
the height range ~ (1.15-1.3) Rq the prominence did not see a PIL with a well 
defined direction. Above 1.3 i?©, the prominence entered the large-scale field 
structure defined by the polar fields and the heliospheric current sheet. The PIL 
direction in this range is relatively similar to the direction of the photospheric 
PIL in the active region and the original orientation of the prominence. Hence, 
throughout the height range of strong initial prominence rotation, ~ (1-1-5) i^©, 
the assumption of rotation by alignment with the relevant PIL is inconsistent 
with the structure of the field. The prominence actually rotated away from the 
PIL in this height range. 

The situation changed when the initial prominence rotation began to level 
off at a large angle. The subsequent rotation up to the height considered in 
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IPatsQurakos and Vourlidasl (j2Qlip is consistent with an alignment with the PIL 
in the hehospheric current sheet. From the orientation reached in the initial 
phase, it was favorable for the prominence to continue the rotation to align with 
the PIL nearly antiparallel to the original orientation. 

These conclusions are robust against the limitations of the PFSS model, which 
reflects the magnetic structure in and around AR 10989 primarily at the time 
of central meridian passage, since the region did not show any strong magnetic 
changes in the time before the eruption on 9 April. The orientation of the PIL in 
the photosphere, as outlined by the prominence in the EUV, changed by about 
20°, much smaller than the observed prominence rotation in the low corona. 
No signs were seen of any major flux emergence which might have altered the 
topology of the ring-shaped PIL seen at low coronal heights. 

Generally, one expects a dominance of the Lorentz force in driving any ro- 
tation of erupting flux at low coronal heights, where the plasma beta is small. 
As summarized in the following section and detailed in Paper II, the Lorentz 
force due to both the tension of twisted field lines and the presence of a shear 
field component in the ambient field causes a rotation away from the PIL, along 
which the rising flux is originally oriented. As larger heights are reached, the 
currents and Lorentz forces in the erupting flux decrease, while the influence 
of the pressure gradient force increases with the increasing plasma beta in the 
ambient field. The alignment with the hehospheric current sheet then reaches 
an increasing importance in the dynamics of the eruption. The transition is 
expected to occur when the plasma beta approaches unity. 

A difference between the PIL in AR 10989 and the large-scale PIL in the 
streamer should be noted. While the PILs are nearly parallel to each other, the 
field direction across them is reversed (Figure [6]). The horizontal field component 
across the PIL in the streamer exerted a Lorentz force on the prominence legs 
according to the mechanism proposed in |Isenberg and Forbes] (|2007p , because 



the prominence was directed approximately along this field component (perpen- 
dicular to the PIL) after the initial rotation had leveled off, i.e., in the whole 
CORl height range. With the current flowing opposite to the axial field of the 
dextral prominence, the direction of this force on the legs of the prominence is 
such that the top part rotates in the clockwise direction. This might explain the 
slight backward rotation indicated by the CORl data. If this were the dominant 
effect in the subsequent evolution, the initial rotation would be reversed until the 
prominence aligns with the hehospheric current sheet in the direction parallel to 
its original orientation. Since the transition to high-beta conditions occurs in the 
C0R1-C0R2 height range, it is unlikely that the Lorentz force could dominate 
the subsequent rotation. Consequently, antiparallel alignment with the PIL in 
the hehospheric current sheet due to the pressure gradient force, which drives 
such alignment on the shortest possible path, independent of the field direction, 
appears to be the most likely evolution in the C0R2 height range. 



3. Numerical Modeling 

Paper II presents a parametric study of various effects which lead to, or influence, 
the rotation of ejected flux ropes in low-beta plasma. The initial flux rope twist. 
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Figure 13. Shape of the erupting flux in the best-matching numerical model — a weakly 
kink-unstable flux rope of initial twist 3.57r — at various stages of the rise, compared with 
STEREO Ahead images at 10:16, 10:26, 10:55, and 11:25 UT. A set of fleld hues enclosing the 
magnetic axis of the simulated flux rope is plotted at the times the rope apex has reached 
the corresponding heights, h = 7.3, 9.6, 20, and 30 ho (ho - initial apex height), which are 
obtained from the scaling of the footpoint distance, 3.3 ho, to the value of 175 Mm estimated 
in Section l2.ll The observed and simulated flux ropes are displayed at the same perspective: 
the line connecting the footpoints of the rope makes an angle of 26° with the line of sight 
and the vertical axis in the simulation is tilted away from the observer by 8° (so that the 
magnetogram, Bz{x,y,0,t), is seen from the bottom side). The dotted line indicates the edge 
of the CORl occulting disk. Since the visible threads in the EUVI images do not outline the 
complete shape of the flux rope, sections of the fleld lines in the corresponding simulation 
snapshots are highlighted for better visual correspondence. For each fleld line, 200 segments 
at the top and 150 segments at each bottom end are plotted with reduced line width. The 
segments are obtained from a numerical integration along the fleld line with adaptive step size 
and are all different, resulting in different lengths of the highlighted sections, analogous to the 
observations. 
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the strength of the external shear field component, and the height profile of the 
overlying potential field are varied. The latter is related to the distance between 
the main polarities in the photosphere. Figure [13] shows field line plots from 
the numerical model that best fits the observations in their entirety (rise and 
rotation profiles and the STEREO images). In this model, a moderate initial 
twist of $ = 3.57r triggers a weak helical kink instability, which lifts the ini- 
tially force- free, toroidal fiux rope into the torus-unstable range of heights and 
also contributes to the rotation. The ejection is primarily driven by the torus 
instability ( [Kliem and Torok, 2006 ), so that the model belongs to the loss-of- 
equilibrium category ( Priest and Forbes, 2002[ ). The major part of the rotation 
in this run is due to the presence of an external shear field component (pointing 
along the prominence and polarity inversion line and due to sources outside the 
fiux rope), as first suggested in Isenberg and Forbes (|2007p . The strength of the 
shear field at the initial fiux rope apex is 2/3 of the external poloidal field which 
holds the rope in equilibrium. The value is consistent with a very rough estimate 
of this ratio from the structure of the active region. The ratio of the distances 
between the main polarities along and across the PIL is about 1/2 in the PFSS 
field, and the shear should slightly increase with the gradual change of the PIL 
direction in the time till the eruption. 

The overall shapes of the fiux rope in the considerable height range included 
in Figure [131 ^is well as the inclination of the still weakly twisted field lines to 
the axis of the rope at the largest height, match the observations reasonably 
well. (Note that the rope has acquired a large part of its total rotation already 
at the first time selected; see Figure [TOl) Similar agreement is demonstrated in 
Paper II for the height-rotation and time- height profiles. 

A model with subcritical initial twist, ^ = 2.57r, which allows only for the 
development of the torus instability, reproduces the rise profile similarly well, 
while the match with the rotation profile is somewhat worse. The field lines at 
the strongly expanded stage of the final CORl image in Figure [13] are rather 
straight and hardly show any indication of the observed twist. Moreover, this 
fiux rope requires a considerable initial perturbation to reach the torus-unstable 
range of heights, which is not supported by the pre-eruption height estimate and 
the initial dynamics of the prominence described in Section [2] (see Paper II). 



4. Summary and Conclusions 

The images of the STEREO EUVI and CORl telescopes resolve the basic struc- 
ture of the erupting prominence in the core of the "Cartwheel CME" on 9 April 
2008 out to heliocentric distances of 4 Rq. A single fiux rope with indications 
of weak to moderate twist is revealed. The fiux rope expands approximately 
self-similarly in the CORl field of view of (1.5-4) Rq. 

The true path of the prominence could be reconstructed from stereoscopic 
observations even when the angular separation between the STEREO spacecraft 
was as much as 48°. The reconstruction reveals that the dextral prominence 
rotated counter-clockwise by a large angle of ^115° up to a heliocentric height 
of 2.5 i?0, where the rotation leveled off. A gentle backward rotation by 15° 
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may have followed in the height range up to 3.3 Rq. Two thirds of the rotation 
were acquired within 0.5 Rq from the photosphere. The coronal height profile of 
the rotation angle in a CME is thus derived for the first time. To our knowledge, 
this ranks as one of the largest rotations so far measured in the corona. 

Taken jointly with recent results about the orientation of the CME at a helio- 
centric distance of ^ 13 Rq, obtained from fitting a croissant-shaped flux rope 
model to stereoscopic images from C0R2 ( Patsourakos and Vourlidas, 2011 ), a 
further counter-clockwise rotation by ^35° from the value at 2.5 Rq is indicated. 
This aligned the erupted flux closely with the heliospheric current sheet above the 
active region. However, the structure of the field within the first ~0.5 Rq above 
the original location of the prominence, obtained from a PFSS extrapolation, 
precludes alignment with the PIL as mechanism for the initial rotation. 

The rotation and rise profiles of the prominence in the EUVI-CORl height 
range are reproduced by a numerical model that follows the evolution of an un- 
stable force- free flux rope. A parametric study of flux rope rotation in this model, 
detailed in Paper II, suggests that the main part of the rotation in this height 
range was caused by the shear field in the source volume. While the rotation 
and rise profiles can be modeled nearly equally well by a weakly kink-unstable 
flux rope of 3.57r initial twist and by a kink-stable flux rope of 2.57r initial twist, 
the indications of twist in the CORl images and the absence of a strong initial 
perturbation in the EWl-Ahead data favor the 3.57r model, suggesting that 
the helical kink instability did occur and contributed the remaining part of the 
rotation. 

In summary, the large total rotation of the erupting flux rope in the Cartwheel 
CME is due to a combination of three effects. The internal Lorentz forces as- 
sociated with the shear field and with the tension of twisted field lines caused 
a strong rotation in the inner corona. The resulting orientation relative to the 
heliospheric current sheet was favorable for the subsequent alignment with the 
sheet in the outer corona and inner solar wind to proceed with the same sense of 
rotation. Since the heliospheric current sheet had an orientation relatively near 
the original direction of the erupting flux, a very large final rotation angle was 
reached. 
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